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Abstract 

We define Sturmian basis functions for the harmonic oscillator and 
investigate whether recent insights into Sturmians for Coulomb-like 
potentials can be extended to this important potential. We also treat 
many body problems such as coupling to a bath of harmonic oscilla- 
tors. Comments on coupled oscillators and time-dependent potentials 
are also made. 

It is argued that the Sturmian method amounts to a non-perturbative 
calculation of the energy levels, but the limitations of the method is 
also pointed out, and the cause of this limitation is found to be related 
to the divergence of the potential. Thus the divergent nature of the 
anharmonic potential leads to the Sturmian method being less acurate 
than in the Coulomb case. We discuss how modified anharmonic oscil- 
lator potentials, which are well behaved at infinity, leads to a rapidly 
converging Sturmian approximation. 

1 Introduction 

A typical situation in quantum theory is to be faced with a physical potential 
V for which one has to solve the corresponding equation of motion, find 
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the energies etc. For most realistic potentials it is impossible to find the 
energies and wave functions analytically, hence one must resort to various 
approximation or numerical schemes. It is the purpose of this paper to extend 
one such approximation scheme from Coulomb-like potentials to harmonic 
and anharmonic oscillators. 

The general situation in quantum theory is a system of N particles interacting 
through some potential V (xi, xn). The equation of motion is then an 
equation of the general form 

(D + V o -E)i> = (1) 

where D is some differential operator of first or second order, Vq is the above- 
mentioned potential and E is an eigenvalue. Examples are (units where 
h = c=l) 

1. The Dirac equation D = i7°7 M c^ (summation over repeated indices 
implied, \i = 0, 1,2,3, 7 M are the Dirac matrices), V = v (x) + 7°m, 
Vo(x) some potential. 

2. The Schrodinger equation D = ^V 2 where V 2 is the Laplace operator 
in d dimensions. 

3. The Klein-Gordon equation where D — □, the d'Alembertian operator, 
V = m 2 + v (x). 

For iV > 1 we simply have D = J2f=i Di where Di is the appropriate operator 
for the i'th particle. Thus all of relativistic as well as non-relativistic quan- 
tum mechanics fall into this category (even in curved spacetimes). Usually 
one finds a complete set of eigenstates corresponding to the different values 
of E, since D is Hermitian these form an ortho normal set. The problem 
is, however, that for many cases continuous as well as discrete eigenvalues 
have to be taken into account in order to get a basis for the full Hilbert 
space. This is for instance the case in the Coulomb potential case already 
in non-relativistic quantum mechanics, where the continuum eigenstates are 
needed to get completeness. Consequently, alternatives will have to be found. 
Shull and Lowdin, [jy , introduced another approach for the Coulomb poten- 
tial. Their methods were generalised by Rotenberg who also coined the word 
"Sturmians" for the new basis set, 0. It has recently been realised that 
this approach can be even further generalised to handle, e.g., many centre 
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potentials, ||] as well as relativistic effects, and many-particle systems [§]. 
It is the purpose of this paper to outline the general theory and to apply it 
to other potentials such as harmonic oscillators and variations thereof. 
Instead of finding a set of eigenfunctions all corresponding to the same "cou- 
pling constants", i.e., charges for the Coulomb and Yukawa potentials, mu 2 
for the harmonic oscillator etc., but to different energies as one normally 
does, one can take the "dual" approach and fix the energy E and then allow 
the coupling constants to vary. Thus we consider not the original equation 
(1) but instead 

(D + p n V -E)i; ri = (2) 

where /3 n is some constant depending on the set n of quantum numbers. 
The solution of this equation gives E as a function of /3 n which can then be 
inverted to find /3 n as a function of E, assuming we can solve the equation 
(0), of course. We will refer to /3 n as the effective coupling constant - for 
Vq = —Z/r it corresponds to scaling the nuclear charge, for Vo = muo 2 x 2 it 
corresponds to scaling the mass and/or the frequency. 
The assumed Hermiticity of D then implies 

((3 n -M fif>* n ,V il> u dx = (3) 

i.e., the Sturmian functions satisfy the potential weighted orthogonality re- 
lation 

/ C^oVveb = N n 5 nn , (4) 

Strictly speaking (Q) only implies / ^,Vo^ n dx = for /3 n ^ /3 n >. For suffi- 
ciently nice potentials, i.e., potentials without any violent oscillations, such 
that the energy E depends monotonically on the quantum numbers n, the 
orthonormality condition ([|) follows. An example of a potential which we do 
not expect to be able to handle with this approach is Vq(x) = 9S ^, but poten- 
tials a; -1 , x~ 1 e~ kx , x 2 can be treated this way. It will also be shown later that 
the many centre analogues of these potentials also are within reach. Another 
subtlety concerns "major" and "minor" quantum numbers in the terminol- 
ogy of Aquilanti and Avery, |5j. The coefficients /3 n need not depend on all 
quantum numbers, those on which it does depend are referred to as "major" 
and the remaining ones are then "minor". For the Coulomb potential, for 
instance, (3 n only depends on n and not on the angular momentum quan- 
tum numbers l,m. The orthogonality with respect to the minor quantum 
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numbers following from the orthonormality of the spherical harmonics Yi m 
and the separation of variables, if)n(r, ft) = Xni(r)Yi m (Q) . For the harmonic 
oscillator no such subtlety occurs. 

To solve the Schrodinger, Klein-Gordon or Dirac equation for some physical 
potential V, we begin by considering Sturmians corresponding to a "base 
potential" Vq for which we can easily solve the corresponding differential 
equation to find the basis set. If Vq is sufficiently similar to V, the conver- 
gence has been found, for the Coulomb and Yukawa potentials, to be very 
rapid, and, if E is taken to be the actual physical energy the Sturmians will, 
furthermore, have the right asymptotic behaviour. We will see, however, 
that this rapid convergence only takes place when the potential (or rather its 
matrix elements) have sufficiently nice convergence properties themselves. 
Equation (||) has far reaching consequences. Consider a new potential V — > 
V = Vq + V', where Vq is some potential for which we can easily find the 
Sturmians (say, Vq ~ r _1 or Vq ~ x 2 ).f\ The equation of motion for this new 
system is then 

(D + Vq + V- E')iP = (5) 
and we can expand ip on the Sturmians for V as ip — En c n^n to obtain 

{D + Vq + V- E') Cn^n = E C " K 1 ~ /W + V + E — E'} ^ n (6) 
n n 

upon using (fj). From this we get the secular equation by using the potential 
weighted orthogonality relation. Thus 

E [(! " Pn)N n 5 nn , + (Wl^Vn) + (E- £')0M^n>] c n = (7) 

n 

from which the kinetic energy has disappeared and only the potential V 
appears together with the overlap integrals of the Sturmians. It is this feature 
of Sturmians which is so important. It often implies that one can find E' 
with very great accuracy from an extremely small basis set, even with just 
one basis function, ip n , in which case we have simply 

E , = (l-/3 u )N n + (il> n \{V , -E!)\i{> n ) 

1 Often, when working with Sturmians, one would like to take V' — V to be a physical 
potential, and only use Vq to find the basis set. We will want to include Vo as part of the 
physical potential in this paper, however. 
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giving an explicit formula for E', which when V is sufficiently close to Vo 
is a surprisingly good fit, [|J [Tj. The reason for this success is to be found 
in the very construction of Sturmians. By construction, Sturmians take the 
potential much more into account, and thus contain much more information 
about the potential. In a sense, Sturmian functions are optimised with re- 
spect to the specific features of the given potential, and it is precisely this 
that lies behind their success in the Coulomb case. 

Another consequence of this can be seen if one attempts to use Sturmians as 
a starting point for a variational calculation. Suppose we know the Sturmi- 
ans for Vo, and now want to use variational theory to estimate the ground 
state energy of the Hamiltonian H = D + Vq + V, using the standard formula 

^ ^ MH\il>) 

En < 



with tj) = if) n , a Sturmian, we get 

which only concerns the starting potential Vo and the energy E to which 
the Sturmians correspond. If the Hamiltonian, moreover, is of the form 
H = D + V', we get 

E <E+ p n TTTT~\ + ii I ; \ W 

which, in any case, is a rather simple integral to compute, suggesting that 
Sturmians be a good starting point for variational calculations. 
It should be noticed, that one will often assume E — E' to be an allowed 
energy also of the "perturbed" potential V = Vq + V for an appropriate value 
of the effective couplings (3 n . In this instance ([!]) can be seen as showing the 
physical effect of the f3 n as a kind of screened charge in the Coulomb and 
Yukawa cases: (1 — /3 n ) measures the "screening". Now, /3 n = 1 corresponds 
to an unsealed coupling constant, from (1 — /? n ) oc (ip n '\V'\ipn) we see that 
for E — E' to be an allowed energy value also for V = Vo, we must fit the 
coupling constant f3 n to exactly the right value as given by the perturbation 
V. 

Since the above consideration only used that D was Hermitian, the results 



5 



hold for all the abovementioned cases provided we can find effective cou- 
plings (5-a such that E can be held fixed. It is known that such (3 n exist for 
the Coulomb, Yukawa potentials both in relativistic (the Dirac equation) and 
non-relativistic quantum mechanics and also for their many centre analogues. 
It will similarly be shown that this also holds for the harmonic oscillator. In 
fact, the following argument seems to suggest that it will always hold. Write 
the original potential as V (x) = av(x) where a is some constant, the original 
coupling constant, specifying the strength of the potential, e.g., a = Z for 
the Coulomb case. Clearly, E = E(a) if E is an allowed energy eigenvalue. 
If there is no degeneracy, a choice for a and for the set of quantum numbers 
n uniquely specify the energy E, and thus we can invert the relation to find 
a as a function of E and n. If there is degeneracy, the equation a = a(E, n) 
has more than one solution in the physical range, but this merely mean that 
more than one set of quantum numbers exists giving the same energy, for 
each such choice of quantum numbers we get a new solution (3 n = a(E, n). 
Another interesting relationship is the momentum space orthonormality re- 
lation. Let (f> n (k) be the Fourier transform of if) n (x), this then satisfies 

(D t - E)(f> n = -/3 n Vo* * n (10) 

where D t , Vq are Fourier transforms of D, V and where * denotes convolu- 
tion, 

V$ * 0n = / V*(k ~ k')<P n {k')dk' (11) 

Multiply by </>*, from the left and perform the /c-integral to arrive at 

J - E)<f) n dk = -(3 n J C,(y * * <j> n )dk (12) 

Now, by the Fourier convolution theorem, (fg) 1 — /* * <?*, and Parzival's 
formula, (f\g) = (f^g 1 ), we get 

J 0*/(-D* - E)<p n dk = -p n N n 5 nn , (13) 

Hence, the momentum space Sturmians satisfy a weighted orthonormality 
relation where the weighting factor is given by the kinetic part of the equation 
of motion and not the potential as in a>space. For the three operators D 
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mentioned in the beginning, this relation reads 



|0 n ,( 7 O 7%-^-7°H0n^ = -f3 n N n 5 nn , 

J </>*, (k 2 - 2mE) <p n dk = -2m/3 n iV n 5 nn , 
J n ,(/c 2 - m 2 - E)(f) n dk = -/3 n iV n <W 

for the Dirac, Schrodinger and Klein-Gordon equation respectively. In the 
non-relativistic case, one will often write k 2 , = —2mE, which is then positive 
(k real) for bound states and negative (k imaginary) for unbound states. 
In this case, the weighting factor becomes k 2 + k% which can be interpreted 
as the length of the momentum vector in d + 1 dimensions. It is this extra 
dimension which is related to Fock's famous treatment of the Hydrogen atom 
where he finds the existence of a S*0 (4)-symmetry. This can be generalised 
to arbitrary dimensions by means of hyperspherical harmonics, []3], Q . 
The momentum space relations have other important implications. The equa- 
tions of motion in momentum space is 

(D t - E)<p n {k) = -(3 n V£ * n (14) 

Define, for simplicity, n = (D l — E)(p n , then we can write the momentum 
space ortho normality relation as / 0*0 n 'C?/c = — /3 n iV n <5 nn / and the equation 
of motion as 

k(k) = -Pn{V3*<i> n ){k) (15) 

= -PnJ V*{k-k')<t> n {k')dk' (16) 

Now, make the following Ansatz 

Vt(k-k')=J2cn4>n(kW n (k') (17) 



Inserting this into (|lq) we then get 



(18) 
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i.e., 



V£(k-k>) = Y.^Mkmk') (19) 

n Pn 

■= E^^X^O)*^)^') (20) 

n Pn 

When Vq is some function of k, k', this is a useful summation formula for the 
momentum space Sturmians. When V * is a differential operator, on the other 
hand, as happens when Vq = x 7 , 7 > 0, then this is a spectral representation 
of that operator. 

It should be emphasised that the coefficients in the above expansion are very 
simple in the sense that they are the natural quantities related to the basis 
set, namely the normalisation, the "effective charge" /3 n , and the kinetic 
operator (which is then a polynomial in momentum space). This shows more 
precisely how Sturmians are adapted to the potential. 
We will now turn to the specific case of a harmonic oscillator. 



2 The Harmonic Oscillator 

To begin with we work with N — 1 and in d — 1 dimension. The potential 
is Vo = \x 2 and we solve the equation 

/ 1 d 2 1 \ 

' +^x 2 - J BU n = (21) 



V 2m dx 2 2 

which is the harmonic oscillator Schrodinger equation with mu 2 replaced by 
f3 n . We will henceforth use mass weighted coordinates and put m — 1, if the 
mass needs to be reinstated one simply replaces x by y/mx. In this case u in 
the original harmonic oscillator Hamiltonian has been replaced by f3^/ 2 - The 
solution of this is clearly 

M*) = 7r- 1 / 4 (n!)- 1 ^2-"/ 2 //„((/3 n )V^)e-^« / ^ 2 (22) 

where n = n is a non-negative integer, H n is a Hermite polynomial and 
E = uj(n + \) = \fb\{n + \). From this relationship between E,n,f3 n we 
read off 

M^) 2 (23) 
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which is the promised relationship between the effective coupling, the quan- 
tum number and the energy, needed to make the Sturmian machinery work. 
The effect of scaling the argument by an n-dependent quantity is to scale all 
the functions to take values within the same interval. This is illustrated in 
figure 1. 

The orthonormality relation, (V'nll^lVv 



^n^nn', then reads 



J HM /4 x)Hni(M 1/4 x)e-^ 1J2 <' 2)x2 x 2 dx = N n ^2 n+1 n\5 nn , (24) 

The first few normalisation factors N n turn out to be as in table 1, we see 
that iV n oc E~ z l 2 - in general, <^ n |a: fc |^n') oc £™( fe+1 )/ 2 for k = 0,1,2,... 
This is easily seen from a simple scaling argument: Perform the scaling 
x — > (3^/ 4 x ~ E 1 / 2 x (for the sake of this argument we can ignore the n 
dependence), in the integral and the result follows directly. 
From table 1 it is clear that 



N = ^ +1 ) 3/2 Jg-3/2 = 1/3-3/4 



(25) 



Notice that the Sturmian orthonormality relation above differs slightly from 
the usual one in two ways, (1) the Hermite polynomials have different argu- 
ments (normally the argument is just -Jujx, irrespective of the value of the 
quantum number n), and (2) the appearance of the factor x 2 . 
Similarly, the overlap matrix between the first five Sturmians T nn / = J ip^ n idx 
is 



/ 



T = E~ 1/2 



i 


1 / 

3 V 



12 /3 
\ 25 V 5 











21 
' 25 



12 

25 



132 
'343' 



21 
25 



fl _132 /30 

343 A/ ""' 







3 

V2 



(26) 



/ 



For d > 1, the energy spectrum is E — u(n + d/2) leading to a very simple 
modification in the expression linking /3 n and E. The wave functions become 
products of Hermite polynomials and n = (ni, rid), n = n i + ^2 + ••• + nd 



A) 



(27) 
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Note, (3 n = (3 n depends only on n = ri\ + ...n^, the "total quantum number". 
The case of more than one particle merely corresponds to a harmonic oscil- 
lator in D = dN dimensions, where N denotes the number of particles and d 
the number of spatial dimensions. If the particles have different masses, one 
have to use mass-weighted coordinates, but otherwise no 

modifications are needed. Clearly, all the important features can be found 
already in the d — N — 1 case, for which reason we will stick to this situation 
in the following unless otherwise stated. 

The momentum space formulation of the harmonic oscillator needs the Fourier 
transform of the Hermite polynomials. It is proven in the appendix that 



^(l + ^n) 1/4 )x 



H n (ik(l + ^ /4 r 1/2 )e^ 2/3 " 1/2 (28) 

is the Fourier transform of ip n . The Hermite polynomials of imaginary argu- 
ments appear in this formula, this means that for n even n is purely real, 
whereas for n odd it is purely imaginary. The momentum space orthogo- 
nality relation then gives the following new relationship between Hermite 
polynomials 

/ H n (ik(l + \^)- 1/2 )H n ,{ik{l + i(/5 n 1/4 )- 1/2 )e^ 2 ^ /2+ ^ 2 )(A; 2 - 2E)dk 
1 n\2 n+1 p 5 J 4 N n (l + '(W (29) 



'2tt 

a somewhat unexpected result 



3 The Anharmonic Oscillator 

We now add a new potential V = ax 3 to V. We then need to compute the 
matrix elements of this in the basis of harmonic oscillator Sturmians. We 
need to compute the matrix W^, = J ifj^x 3 ifj n >dx. For n, n' = 0, 1, ...,4 we 
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get the following explicit result 



= E~ 2 



V 



o 

27 
64 



343^ 
2048 





27 

64 



7425^ 
2048 



3159 
512V 7 ! 





7425^5 
2048 



677425^ 
41472 





343 \/3 
2048 



677425^ 
41472 



431831169 
4194304 





3159 
512^ 



431831169 
4194304 





(30) 



The secular equation then reads 



det((l - /3„)iV n (W + aW (3) ) = 



(31) 



By including the first N Sturmians we get the ground state energies shown 
in table 2 with a — .1. 



Another version is V 
becomes 

/ 



x A , for this case the matrix W, 



(4) 



W (4) = E -5/2 





25 
144 



243 
1000 



^ 



135 
16 



27783 
2000 





25 
144 



975 
16 



625725 
9604 





27783 
2000 



3675 
16 





243 
1000 



625725 
9604 



29889 



3 \ 



/ 



1. 



16a/2 

(32) 

, 5 Sturmians 



As for x 3 we get the energies E by including the first N 
also shown in table 2. As is apparent from the table, the ground state 
values are unstable, i.e., the approximation has failed to converge. This is 
in sharp contrast to what is known to be the case for the Coulomb and 
Yukawa potentials (both in their single as in their many-center form) where 
the convergence is very rapid. This failure can be traced back to the non- 
convergence of the matrix elements W^, — > oo for n,n' — > oo, which again 
is a consequence of the divergent behaviour of the potentials as x — > oo. 
Although the Coulomb and Yukawa potentials are singular at the origin, r = 
0, their matrix elements none the less fall of rather rapidly as the quantum 
numbers increase, due to the quick fall-off of the potentials themselves as r 
increases. 

Standard perturbation theory would give a value for the ground state energy 
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in the two cases of 



11 h 2 a 2 

EL = E 5— (cubic potential) (33) 

8 m 6 ijj^ 

. 3 h 2 a 23 ha 2 

E = E + — — - (quartic potential) (34) 

16 rrrur 4 m A uj 6 

by going to second order in the coupling constant e and using standard 
Rayleigh-Schrodinger perturbation theory, which gives a non- convergent se- 
ries - for the particular example of m = h = u = 1, a = .1 we get E' = 0.3625 
for the cubic and E' Q = 0.46125 for the quartic anharmonic oscillator. Clearly, 
this is not in good agreement with the result found by using Sturmians, but 
as will be seen below, this is due to the divergence of the perturbation se- 
ries. That the Sturmian method is non-perturbative is suggested by the 
general solution of the secular equation, for the quartic anharmonic poten- 
tial the ground state energy as a function of a is found to be (for N = 1) 
E' = z 1/3 + jzZ- 1/3 where 



2592a + v / 6718464a 2 - 6912 

Such a dependence of E on a can only be obtained in perturbation theory 
by performing at least a partial resummation of the infinite series. That 
the perturbation series is divergent is mirrored in the behaviour of E' as a 
function of the coupling constant a: As a increases so does the real part of 
E' Q . Furthermore, E' has a small imaginary part which is wildly oscillating 
but vanishes for a > 0.003 and hence quickly becomes unimportant. 
Thus, we must compare our Sturmian energies with non-perturbative results. 
The divergence of the perturbation series can be improved in a number of 
ways, as done recently by Bender and Bettencourt, ||, and by Kunihiro, 
||. The latter paper performs a re-summation of the perturbation series 
by using a renormalisation group (RG) improved technique. From this very 
rapid expressions for the ground state energy for the quartic anharmonic 
oscillator is found. Some high-precision numerical results for the quartic 



anharmonic oscillator have been found by Bacus et al., |TIJ, and we will make 



a comparison with those findings. Other recent papers on related topics are 



11 



We will go back to the original secular equation, (0), and let E be the energy 
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of the harmonic oscillator. When E = E n = (n + l/2)Hu we get /3 n = 1 
(we will say we are "on shell") and the secular equation simplifies, since the 
results then no- longer depend on the normalisation factor N n oc E^l 2 . Let 
the number of Sturmians in our basis set be N, then for iV = 1 we get simply 

E' = E + JJ^M = E + cO£w!& (35) 

Let N = 2 and let the basis set correspond to the quantum numbers n, n + 1, 
we then find E' by solving the quadratic equation 

= (E-E , ) 2 T nn T n+hn+1 + a(E-E')(T nn wi% tn+1 + 

T n+1 , n+1 W^) + a 2 W^wi% n+1 (36) 

For n = 0,1,.., 4 we get the results shown in the table 3, where we also 
show the highly accurate (approx. 90 significant digits) of Bacus et al.. We 
note that this crude Sturmian approximation is able to get the right order 
of magnitude for the energy states even of the excited states, but we also 
notice that the accuracy decreases as n increases. This is once more due 
to the non-convergence of the matrix elements of the pertubation potential. 
Increasing the basis set will not lead to improved accuracy, but can in fact 
lead to quite the opposite because of this divergence. Hence, the best results 
are found by the simplest approximation, namely N — 1. In contrast to 
the case of asymptotically well-behaved potentials such as the Coulomb and 
Yukawa ones, hitherto studied in the literature. 



4 A Comment on Time Dependent Poten- 
tials: The Damped Oscillator 

Another important variant of the harmonic oscillator is the damped os- 
cillator. The particular Hamiltonian we're going to study is the so-called 



Caldirola-Kanai oscillator, 12 



H = V -e~ 2 ^ + -ux 2 e 2 ^ (37) 
2 2 
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where 7 is some constant, the friction the coefficient and t denotes time. This 
can be re-expressed in terms of a potential V 

V = -cuVe^'W = Z(t)V (x) (38) 

with 

1 - e~ 27 * 

fit) = (39) 

7 

For time dependent potentials V, the secular equation have to be modified. 
The Schrodinger equation for the full system reads 

(d + v + v'(t)W = ih^ (40) 

at 

expanding ip = J2n c n(t)ip n leads to the modified secular equation 

E K 1 - Pn)N n 5 nn , + (MV\t)\il> n ) + E(^ n }} c„ = i^c n (i-|i) 

n n 

(41) 

Using the relationship between V and V we can rewrite this as 

E K 1 "^(1 +eW)Wnn' +£(VV^n>]c„ = ^^^{^'1^) (42) 
n n 

For the simplest possible case of only including N — 1 Sturmians in the basis 
set, the solution to this secular equation is of course 



/ r A/ 1 TV /■* \ 

c n (0=c n (0)exp(-i* + E -Hi-fo-JE- ^ t ')dt') (43) 

V L inn J inn ^0 / 

It turns out that one can actually compute the integral, since 

and the integral can be expressed in terms of the exponential integral function 
leading to 

c n (t) = c n (0)e- iuJot - iuJ ^ (44) 



where 



N 

uio = (l-/3n)-± + E (45) 

J- n 



nn 

o2 



u(t) = (l-/3 n )^|-(Ei(-2)-Ei(-2e- 2 ^) (46) 

i nn ^ I 
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Naturally, "on shell" (i.e., for /3 n = 1) we simply get c n (t) = c n (0)e~ tEt as 
one would expect, whereas "off shell" (i.e., for (3 n ^ 1) we get a highly 
oscillatory behaviour. The explicit results for iV = 1, n = are u = 
±±g^,u;(t) = e2{ {~^ 2) (Ei(-2) - Ei(-2e" 2 ^)). One should also notice that 
this result holds even for N ^ 1, one must then interpret the division by T nn 
as multiplication from the left by the inverse T~^,, whereby u>o,u(t) become 
matrices. 

The real and imaginary parts of Co(t)ipo(x) have been plotted in figure 2a-b. 
We notice that the oscillations in the t-direction die out either as x increases 
(because of the decrease of ipo(x)) or as time goes. 

This simple example shows how time-dependent problems simplify in the 
Sturmian approach because of the simplification of the secular equation. 
Thus, Sturmians are well suited for problems with time- dependent poten- 
tials or for scattering processes. But they are of course subject to the same 
limitations as in the time-independent case. 

It also shows that this particular time-dependent damped oscillator is exactly 
solvable using the Sturmian approach. 



5 A Bath of Harmonic Oscillators 

Consider the potential 

V = W^~^) 2 (47) 
z i 

where i runs over some index set. This represents the potential coming from a 
family, indexed by i, of harmonic oscillators situated at Xj and with coupling 
constants (characteristic frequencies) g^. We will usually restrict ourselves to 
i being discrete corresponding to an infinite lattice of oscillators, but for field 
theoretical purposes it can also be relevant to allow i to run over a continuous 
index set (in which case the sum must be interpreted as an integral). This 
example is the harmonic oscillator analogue of the many centre Coulomb 
potential treated in j|, ||. Note that we can rewrite V as 

V = ^gx 2 -xJ2 9i^i + - 9&i ( 48 ) 

i i 

= gVo-x^ + 02 (49) 
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where g = J2i9i, c i — Hi9i x ii c 2 = ^J2i9i x h hence we can see this as a 
perturbation of the original potential Vq. Because of this feature of the 
bath of harmonic oscillators, our computations will simplify somewhat from 
the many-centre Sturmians introduced for the Coulomb potential in ||, [|], 
which is very fortunate since those papers use a Fourier transform approach 
which is not useful for the harmonic oscillator - as mentioned earlier, the 
Fourier transform, V t (k), of the potential V is a differential operator for 
the harmonic oscillator (actually, the Hamiltonian is form-invariant under 
Fourier transforms). 

Furthermore, we can complete the squares to write V as another harmonic 
oscillator plus a constant, in fact 

„ \ 2 „2 

(50) 

Consequently, we can obtain a solution to the many-centre Schrodinger equa- 
tion by simply making the following substitutions in the solution for the single 
harmonic oscillator 

2 

x^x + ^ P n ^gP n E^E-c 2 + ^ (51) 




i.e., the new /3 n , (3 n reads 



Pn = 9\ 



71+1/2 



Hence the Sturmians read 



* n (x; x t ) = TT-V^n!)- 1 /^-^ Ifc/* {x + | } j e -Wh(*+%f (52) 

These are then a convenient basis for many-centre problems. 
For the many-centre Coulomb potential, the Sturmians become related to 
the one-centre Sturmians by means of phase factors e tkiX \ whereas for the 
harmonic oscillators the many-centre Sturmians are related to the one-centre 

c 2 

ones by means of a translation well as a scaling and a shift in 

energy. 

As a simple example, consider a particle moving in a potential coming from 



16 



harmonic oscillators situated at Xi = i,i = 1, ...,M all with equal strengths, 
& = I. Then g = M, a = §M(M + 1), c 2 = ^M(M + 1)(2M + 1) and 

x^x + ^(M+l) E -> E - ^M(5M 2 + 6M + 1) 

thus 

~ /^ 96E - 2M(5M 2 + 6M + 1) \ 2 

^ n " ^ 48(2n + l) J 

and the first two Sturmians read explicitly 

t/ n i/4 { »A8E-M(5M 2 + QM + l){ \. %r A 2 \ 
Voix-Xi) = Tr-^expf-M ^ ^ (x+-(M + l)j52j) 

= tt- 1/4 ^ v /m(48E - M(5M 2 + 6M + 1)) + -(M + 1)) x 



exp f -M K — l - [x+-(M + 1)J(&) 

The energy is found from the Schrodinger equation which leads to the "on 
shell" condition f3 n = 1. Consequently, the energy of the n'th state is 

E n = (n + ^)F 1/2 - c 2 + ^ (55) 

irrespective of the number of Sturmians used as the secular equation (by 
construction) is diagonal. In the particularly simple case of M evenly spaced 
oscillators all with the same value of the coupling = 1 readsQ 

En = (n + ^M- 1 / 2 - —M(M + 1)(2M + 1) + j-M(M + l) 2 (56) 
2 12 lo 

which is then the energy of a (non relativistic) particle moving in a one- 
dimensional lattice of oscillators - a highly simplified model of, say, a particle 
in a solid. For a many- dimensional lattice we would simply use products of 
unidiemensional Sturmians. 



2 Note, for M = 1, E ^ n + 1/2, because the potential in this instance is V = \{x~ iy 
and not V = Vq. Had we instead used Xi = i — 1, we would get V(M — 1) = Vq, E(M = 
1) = n + 1/2. In that case, by the way, c x = \M(M - 1), c 2 = j^M(M - 1)(2M - 1). 
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6 A Comment on Coupled Oscillators 



Consider now a potential of the form 

V=^ y Z,9i(x-Xi) 2 + \T, X ij( x - x i) 2 ( x - x j) 2 ( 57 ) 

i i^j 

which introduces a coupling between the oscillators at the various positions. 
In a manner similar to the manipulations of the bath of oscillators, this can 
be transformed into a single anharmonic oscillator potential 

V = -gx 2 — c\x + c 2 — c 3 x 3 + c 4 a; 4 (58) 

where 

9 = Eft ( 59 ) 

i 

ci = E 2 (9i - E + x i x i)) x i (60) 

c 2 = ^E(ft + E^ 2 K 2 (6i) 

i j^i 

c 3 = ^2Xij(xi + Xj) (62) 

c 4 = \Y. X H (63) 

Suppose x is a solution to the corresponding Schrodinger equation. We can 
then expand x either on the ordinary Sturmians, ijj n , or the Sturmians for a 
bath of oscillators, \I/ n . If we choose the latter option, we have to complete 
the squares to obtain the centre of the new oscillator, but this would mean 
that the anharmonic terms too would have to be shifted and this would 
again introduce lower powers of the new, shifted position. Hence the secular 
equation would end up having the same structure and thus the same level 
of complication. Consequently, nothing is lost by expanding on the single 
oscillator Sturmians tp n , x — En Q n^n^ The secular equation then becomes 

= E [(& + g)NJ nn , - Cl w£>, + c 2 T nn , - c 3 W® - c 4 W^] cv (64) 
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where 



W, 



ip^x k %jj n idx 



(65) 



is the matrix elements of the fc'th power of x, W^, = T nn i. The only one of 
these we do not already know is for k = 1 in which instance a straightforward 
computation yields 



/ 



= E- 1 



\ 





3 
8 



21 /o~ 
128 V" 3 





3 
8 



75 /o~ 
128 V<3 



81 
64^ 





|v 




128 ^ 



1925 /F 
1728 v 





_ 21 /q 

128 V" 3 



1925 /F 
1728 v 



508599 
131072 



\ 



81 



64^2 


508599 
131072 





(66) 



/ 



for the first five Sturmians. Consequently, for N = 1 we get (since W r< - 2fc+1 - ) 
is diagonal) 

(Pn + 9)N n + c 2 T nn -c 4 W^ = (67) 
which is a cubic equation for E, e.g. for n = 



g - 3c 4 E + AE 2 + 16c 2 E 3 







which for the extremely simple case of g = c 2 = C4 = 1 has the three solutions 
E = -0.669498, 0.209749 ±i0.222168, i.e., one negative energy state (hence 
a bound state) and two complex conjugate oscillatory states. 
For N = 2 and g = C{ = lQleads to the following solutions E = —9.91107, —1.51155 
and E = 0.129506 ± iO.435961, 0.537995 ± z'1.32394, which then corresponds 
to two bound states and two pairs of oscillatory states, the latter of which 
essentially oscillates around the ground state of the single harmonic oscilla- 
tor. 



7 The Gaussian-Damped Anharmonic Oscil- 
lator 

The previous computations seem to suggest that the Sturmian method is best 
suited for potentials which are well behaved at infinity such as the Coulomb 

3 Which, by the way, is only possible for two coupled oscillators if x\ = 4(1 =F 
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potential, but very slow converging for potential diverging as x — > oo such 
as the anharmonic oscillator. To test this hypothesis, we will now briefly 
consider a toy model, the anharmonic oscillator damped by a Gaussian V = 
x where k — 3, 4. 4 Including only the first five Sturmians we get the 



ax k e 



following matrix elements, W^> of V, k = 3, 4 



W {3) 



E 1 ' 2 





27 



27 



2(3+4E)5/ 2 





3 343(3+2E) 
2 2(7+8£) 7/2 







2(3+4_B)5/ 2 



225(15-22£) 

4 v / l+¥( 15 + 8E ) 3 



81(243-324g-52g 2 ) 
4 v /2+§S(9+4E) 4 





225(15-22E) 



3 343(3+2£) 
2 2(7+8£) 7 / 2 





1225(525-94Qg+316g 2 )>/21 
4 v /l+|fB(35+12 J B) 4 



250047(107163-207522E+94356E2_i38l6E 3 ) 
8(63+16£) n /2 



4 v / F¥( 15 + 8S ) 3 



1225(525-940E+316-E 2 ) V ^T 
4 % /l+i§B(35+12E) 4 







4(l+2£) 5 /2 



75(5-4E) 





243(243- 160E 2 ) 

V 16 v / l+¥( 9+10E ) 4 





5 

2(l+2E/3) 7 /2 



19845g(21-4g) 



^£(21+10^)4 







81(243-324E-52g 2 ) 
4 v /2+|E(9+4£) 4 



250047(107163-207522£;+94356£;2_i38l6 J E 3 ) 
8(63+16B) 11 /2 







75(5-4£) 



(68) 



4 One could also consider simply the exponentially damped 
due to the presence of the absolute value the matrix elements 





75(25-8Qg+104E 2 ) 
8 A /l+§E(5+4S) 4 

2025(273375-899100-E+1219680g2-92288£; 3 ) 

16 v / FW( 45 + 14 ^) 5 

oscillator V = x k e~\ x \, but 
become too complicated. 
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243(243- 160-E 2 ) 



\ 







19845£(21-4g) 



16^1 + ^(9+10^)4 





2025(273375-899100E+1219680g 2 -92288E 3 ) 



(69) 



735g(147-112g+40E 2 ) 
4Sgrtl+^(7+2£) 5 








/ 



which leads to the ground state energies shown in table 4 below. We notice 
the improved convergence properties supporting our claim that it was the 
non-convergence of the matrix elements of the undamped anharmonic oscil- 
lator that was the cause for the non- convergence of the ground state energies 
in the Sturmian method. 

8 Conclusion 

We have seen that the powerful technique of Sturmians functions developed 
for Coulomb-like potentials can be extended to harmonic and anharmonic 
oscillators, where it furthermore can be seen that the technique is highly 
non-perturbative, but the divergence of the potential (as x — > ±oo) leads to a 
non-convergence of the Sturmian approximation, in contrast to the Coulomb 
case, where we have very rapid convergence. It turned out, however, that 
already with N — 1, 2 Sturmians the correct order of magnitude for the en- 
ergies of even the excited states could be obtained. Thus indicating that the 
problem with convergence is perhaps not so serious after all, if one merely 
wants to find the order of magnitude. For higher precision, one should prob- 
ably utilise a hybrid method, using the first few Sturmians to get the correct 
order of magnitude and then some variational approach, say, to get the re- 
quired precision. 

We also saw how to treat time-dependent problems, where once again the 
Sturmian properties lead to some important simplifications. Finally, we con- 
sidered a bath of coupled or uncoupled oscillators which could be transformed 
into a single anharmonic oscillator problem. This is contrary to what one 
does for the Coulomb potential, where Fourier transform techniques are used 
in stead. Furthermore, when using a modified (or regularised) potential, con- 
vergent at infinity but with the same behaviour for x not too large, we did 
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get rapid convergence, especially for the x 3 case, whereas the x 4 case had 
slightly slower convergence. 

All of this seems to suggest that the Sturmian techniques have a very wide 
range of applicability covering basically all important potentials known in 
atomic physics or quantum chemistry, but one should be very careful when 
using potentials which are not well behaved at infinity. Given the general- 
ity of the approach, as outlined in the introduction, this procedure should 
also be extendible to problems in quantum field theory using a functional 
Schrodinger picture and to problems in quantum kinetic theory in phase 
space using Wigner functions for instance. 
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A Certain Results Concerning Hermite Poly- 
nomials 



In this appendix we prove a few results concerning Hermite polynomials. The 
generating function is known to be 



from this it is straightforward to deduce the standard orthonormality relation 
for the Hermite polynomials and the harmonic oscillator wave functions. We 
will also use it to derive some other useful results. 

First we need the Fourier transform of a harmonic oscillator Sturmian. Thus 
we want to compute the following integral 



The corresponding integral with the generating function is merely a Gaussian 
integral and can be readily computed 



oo n 



e ~s 2 +sx _ 



£ -,H n {x) 

n=0 ,L " 



(70) 




st+asx-lpxi+ikxfa = e -s 2 f^ e jp(as+ik) 2 




(71) 
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= £ -7 / H n (ax)e-^ x2+lkx dx (72) 

n=0 n - 

from which we get, by Taylor expansion, the desired result: 

J H n (ax)e-^ 2+ ^dx = JjU + ^(^(iW/f^))- 1 ^^ 1 ^ 

(73) 

Actually, this formula is a little more general than we need. For harmonic 
oscillator Sturmians it turns out that a — (5 (— {P^m) 1 ^) which leads to a 
slight simplification, resulting in the formula given in the text. 
Another expression we need is the matrix element of x 7 for 7 some positive 
integer > 2, i.e., we need to compute 

I nm (a,P,>y,6):= J ' H n {ax)H m {(5x)x^ 'e~^ 'dx (74) 

Again, in terms of the generating function, the integral we need to compute 
is quite simply 

/V 2 « a+2 <° + ' 3 >"-i fc V<ir = e-^V^r'-^ra + 1) x 

(2<„ + /jMi-(-in lFl(1+ i.i2^±f^ )+ 

(75) 

= l(a,P,y,6) (76) 

which is valid even for non-integer 7. Here \Fi is a hypergeometric function 
and L b a is an associated Laguerre polynomial. The first few of these are 

Li 1/2 {x) = ^-x 

L 2 1/2 (x) = -(3-12x + 4x 2 ) 
8 

Lz 1,2 (x) = ^(15-90a; + 60a; 2 -8x 3 ) 
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i^4;^) = ^(3 + 2x) 

= ^e*(15 + 20x + 4x 2 ) 

For the two cases of interest to us, 7 an even or odd positive integer, we get a 
simpler relation since one of the two terms on the right hand side will vanish. 
For 7 = 2k we arrive at 



l(a,(3,2k+l,5) = 2 k ^e^ 2 5- k -^T(k+ 3 -)s(a+P) 1 F 1 (k+ 3 -; §;2 (a + .f )V 

2 2 2 



whereas we for 7 = 2/c + 1 find 

3 :"; in D 2 .k- 

Taylor expanding in s, we get 

00 00 ^n+m 

n=0 n! n,m=0 n!m! 

From this we can obtain relationships between the matrix elements W^> f° r 
different fc's, since 7 nm oc for a = = fiH^- For instance 

-^00 — Z 

I01 + -Tlo — 2Ti 
f 1 1 
2^02 + 2^20 + ^11 = t^ 2 

^(/03 + /30) + ^(/21+/l 2 ) = 

such rules (essentially following from the recursion relation for the Hermite 
polynomials) can be used to simplify the computation of matrix elements. 
One should also note that, since the above is valid even for 7 not a positive 
integer, we can use it to obtain the matrix elements of the Coulomb potential 
between harmonic oscillator Sturmians. 
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n 


N n 


n 





1 

2 
3 
4 


_J_ £-3/2 

9/3 £-3/2 
8 V 2 

25 fbp-3/2 
8 V 2 

49 ^/Z 77-3/2 
8 V 2 
243 R -3/2 


5 
6 
7 
8 
9 


±§y^- 3/2 

169 /Up-3/2 

8 V 2 
225 /Up-3/2 

8 V 2 
289 /17 p-3/2 

8 V 2 
361 fWp-3/2 

8 V 2 



Table 1: Normalisation constants for harmonic oscillator Sturmians. 



N 


£(:r 3 ) 


E(x 4 ) 


1 


0.500000 


0.562709 


2 


0.014628 


0.562709 


3 


0.112767 


0.562709 


4 


0.351135 


0.562544 


5 


0.102981 


0.562516 


10 


1.27012 


0.533858 



Table 2: The energies for the cubic and quartic anharmonic oscillator V = 
with a = .1 found by using only the first N Sturmians. Only the 
ground state energies are shown. 
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n 


E (Sturmian) 


E (Bacus) 


difference 


N = 1 


N = 2 


N = 1 


N = 2 





1.07500 


1.07500 


1.06529 


-0.00971 


-0.00971 


1 


3.37500 


3.37500 


3.30687 


-0.06813 


-0.06813 


2 


5.97500 


5.97500 


5.74795 


-0.22705 


-0.22705 


3 


8.87500 


7.00152 


8.35268 


-0.52232 


1.35116 


4 


12.0750 


9.30093 


11.09860 


-0.97640 


1.79767 



Table 3: A comparison between the perturbed energy states for the quartic 
anharmonic oscillator with a = .1 found by using two Sturmians and the 
high-precision results of Bacus et al., the second order perturbative result for 
this case (m = \,u = 2) is E = 0.4900 for the ground state. 



N 


E(x A ) 


E(x^) 


1 


0.500000 


0.622877 


2 


0.495852 


0.622877 


3 


0.491822 


0.622878 


4 


0.491282 


0.622877 


5 


0.491282 


0.622878 



Table 4: The ground state energies for V = x k e x , k — 3, 4 computed using 
N Sturmians. 

Figure 1: A plot of the first ten harmonic oscillator Sturmians showing how 
they are scaled to all have essentially the same range. 
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Figure 2: The real (a) and imaginary (b) parts of the first time depen- 
dent Sturmian for the damped (Caldirola-Kanai) oscillator, c n (t)ipo(x), in 
the range t G [0, 10], x e [0, 5] and with E = 7 = 1. 
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